SUNY-NTG-95/50 



Scalar susceptibility in QCD and the multiflavor Schwinger model 



A. Smilga 

ITEP, B. Cheremushkinskaya 25, Moscow 117259, Russia 

and 

J.J.M. Verbaarschot 

Department of Physics 
SUNY, Stony Brook, New York 11794 

Abstract 

We evaluate the leading infrared behavior of the scalar susceptibility in QCD and in 
the multiflavor Schwinger model for small non-zero quark mass m and/or small nonzero 
temperature as well as the scalar susceptibility for the finite volume QCD partition func- 
tion. In QCD, it is determined by one-loop chiral perturbation theory, with the result that 
the leading infrared singularity behaves as ~ logm at zero temperature and as ~ T j ^fm 
at finite temperature. In the Schwinger model with several flavors we use exact results for 
the scalar correlation function. We flnd that the Schwinger model has a phase transition 
at T = with critical exponents that satisfy the standard scaling relations. The singular 
behavior of this model depends on the number of flavors with a scalar susceptibility that 
behaves as ~ j72~2/(a'/+i)_ flnite volume V we show that the scalar susceptibility is pro- 
portional to Xjm^V . Recent lattice calculations of this quantity by Karsch and Laermann 
are discussed. 
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1 Introduction. 



The scalar susceptibility in QCD is defined as 

Nf Nf Nf 

X= d'x{Y.(lMx) Y^q.q^iQ)) - V{Y.m^? = y dl\ogZ\^^^, (1.1) 

i=l i=l 2=1 

where V is the four dimensional Euclidean volume and the averaging is performed either 
over the vacuum state or over the thermal ensemble. (In the Euclidean approach the 
latter corresponds to an asymmetric box with imaginary time extension of /5 = 1/T ^ L, 
i.e. V = L^P). The definition (|1.1| ) is for a diagonal mass matrix with equal quark 



masses. It is especially interesting to study this quantity in the neighborhood of the 
thermal phase transition point. It is expected that for QCD with two massless flavors 
a second order phase transition occurs leading to restoration of chiral symmetry [|| (see 
0, ^ for recent reviews) with a diverging susceptibility. This can be understood simply 
in terms of Landau mean field theory. For a system with order parameter t] coupled to 
external field h with a second order phase transition at T = Tc, the fluctuations of the 
order parameter are described by the effective potential 

V^^{7]) = A{T-T^)r]^ + B7]'^ + C7]h. (1.2) 

In QCD 1] is the chiral condensate and h is the quark mass. At T = Tc, the minimum 
of the potential occurs at ~ h^^^ which gives the law x = drj/dh ~ for the 

susceptibility. On the other hand if T is close to but T ^ and the external field is 
weak enough 

/K|r-T//^ (1.3) 

the quartic term in (|1.2| ) is irrelevant and the scaling law is x ~ \T — Tc\^'^ both above 

and below Tc (the proportionality constant in these two regions differs by a factor 2). 

Recently, on the basis of lattice simulations of the 3-dimensional Gross-Neveu model 

it has been suggested that a second order phase transition involving soft modes consisting 

out of fermions has critical exponents given by mean field theory . In particular, if this 

is also true for QCD with two massless flavors, we get 

1 

~ m3, (1.4) 

X ~ (1.5) 



at T = Tc, and 



(1.6) 

at 

A^mi < |r - Tel < T„ (1.7) 

where A is a typical hadronic mass scale. These scaling laws have been reproduced by a 
simple stochastic matrix model 0, The scalar susceptibility was recently measured for 
lattice QCD with two light flavors |^. They found a diverging susceptibility at T = with 
critical exponent = 0.24 ±0.03 which does not agree with the prediction = 1/3 of 
the Landau mean field theory. Obviously, further numerical measurements of the critical 
indices in QCD are highly desirable. 

In this paper we will study the scalar susceptibility in three different cases. First, 
in section 2 the scalar susceptibility is evaluated for the multi-flavor massive Schwinger 
model which shares many common qualitative features with QCD and, in particular, it 
shows a phase transition with "restoration" of chiral symmetry at T = 0. All other critical 
exponents that can be defined in the SM will be evaluated as well, and it will be shown 
that they satisfy the scaling relations modified for a phase transition at zero temperature. 
Second, we evaluate the scalar susceptibility in QCD at low temperatures using chiral 
perturbation theory. Third, in order to estimate finite size effects in lattice calculations, 
we calculate the scalar susceptibility in volumes with spatial length below the Compton 
wave length of the pion. 




a b 



Fig. 1. Connected (a) and disconnected (b) graphs contributing to the scalar susceptibility 
of quarks propagating in a background gluon field. 

Before proceeding further, we should note that the scalar susceptibility as defined 
in (|l.lj ) involves a quadratic ultraviolet divergence due to a trivial perturbative graph 

3 



depicted in Fig. la. The situation is the same as with the chiral condensate which 



involves a trivial divergent perturbative contribution ~ mAfjy. The predictions (|1.4|JT75 



L6|) hold for the infrared sensitive part of {qq) and x- 

The susceptibility (|1.1|) can be written as the sum of a connected and a disconnected 
contribution corresponding to the graphs of Figs, la and lb, respectively, 

X = Nfx''^^ + N]x^'^ (1.8) 

The connected contribution to the susceptibility was calculated in The disconnected 
contribution is defined by 

= J d^x{uu{x)dd{Q)) - V{uu){dd) = ^dm^dm, \og Z , (1.9) 

where and are two different quark masses that are put to zero after differentiation. 
We will obtain the latter contribution from the difference of x and x'^"'^- 

Both scalar susceptibilities can be expressed in terms of the eigenvalues of the Dirac 
operator 

V 2\ / ^ V 2' 



ydis ^ 1 
A y 



;i.io) 



A related susceptibility is the pseudo-scalar susceptibility which in terms of the Dirac 
eigenvalues is given by 



2 The Schwinger model 

In this section we discuss the Schwinger model (SM) with Nf light flavors, and determine 
the critical exponents from known results for the correlation functions. For the critical 
exponents we will follow the convention of f^. 
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Critical behavior shows up only in the SM with several flavors. In the standard 
Schwinger model {Nf = 1), there is no non-anomalous global symmetry to be broken 
spontaneously and no reason for the phase transition to occur. The scalar susceptibility 
in the standard SM was determined recently in [|l^. It is just a finite constant. 



The theory with Nf > 1 and zero fermion masses has the global symmetry SUL{Nf) x 
SUji{Nf) (much like in QCD) and the potential possibility of its spontaneous breaking 



with generation of fermion condensate exists. Coleman's theorem |]TT| prevents, however, 
spontaneous breaking of a continuous symmetry in 2 dimensions. So, a QCD-like phase 
transition cannot occur in a 2D theory at a nonzero temperature. Nevertheless, it turns 
out that the dynamics of the SM at small T is similar to that of a theory with second 
order phase transition in the region of temperatures slightly above critical. One can say 
that the phase transition does occur at T = 0. 
The Lagrangian of the model is 

^ = -jPlu + ^ ^/7m(^m - i9K)qf -mY, QfQf, (2-1) 
^ /=i /=i 

where g is the coupling constant with the dimension of mass, and, for simplicity, we 

assumed that all quark masses are equal. 

Let us study this model in the region m <^ g. The particle spectrum involves a massive 

photon [jl2| with mass 

/i+~(?J^ + 0(m) (2.2) 

V 71 



and 'quasi-Goldstone' particlesQ with the mass [T^, [14 



^_ ^ (2.3) 

This gives us the critical exponent fi = Nf/{Nf + 1). 

At nonzero m, the chiral SUiiNf) x SUnlNf) symmetry of the massless SM Lagrangian 
is broken explicitly, and the formation of the chiral condensate becomes possible. The 
chiral condensate involves a UV divergent piece 

-mlogA[/y, (2.4) 



"'^It is better to use the term 'quasi-Goldstone' than the term 'pseudo-Goldstone' commonly used for 
pions because 'quasi-Goldstone' states in the SM become sterile in the chiral limit ni — > 0. That conforms 
with Coleman's theorem which forbids the existence of massless interacting particles in 2 dimensions. 
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and an infrared contribution (sensitive to the small eigenvalues of the Euclidean Dirac 



operator). The latter has been determined in 0, ^ with the result 



{qq) ^m''f+'g''f+\ (2.5) 
providing us the critical exponent 6 = [Nf + 1)/(A^/ — 1). The susceptibility is 



X 



(-) (2.6) 



dm 

In the region /i^^ ^ |a;| ^ fiZ^ the vacuum scalar correlator is given by []T2|, |T^, |T^ 

(gg(x)gg(0))o ~ jS^^- (2.7) 



\x 



The associated critical exponent is C = 2 — 2/Nf (To make contact with the standard 
theory of phase transitions, x should be assumed space-like, but the behavior (|2.5| ) holds 
of course for any x due to Lorentz-invariance ). At ^ fiZ^ the correlator levels off at 
the value of the square of the chiral condensate ( p.5|) . 



In the region ji^ <^T <^ ji^^ (weak field limit) the condensate is given by ||T5 

'9 



2/Nf 

(qq) ~ m ( I ) , (2.8) 



which leads to the susceptibility 



X ~ j , (2.9) 



and the critical exponent 7 = 

The correlation length in this case is just inverse fermion Matsubara frequency ~ 1/T 
which gives the critical exponent u = 1. For <^T the energy density is that of a gas 
of massless particles. Therefore, e ~ and the specific heat for zero field is given by 

C = ^^T (2.10) 

leading to a = — 1. 

The above results for the critical exponent have been summarized in Table 1. We also 
show the critical exponents for mean field theory. Note that for Nj = 2 some (namely, S 
and 7), but not all of the exponents coincide with the predictions of MFT. 



exponent MFT SM 
a ~ 



7 1 
5 3 



2 

Nf 
Nf+1 



— 

1 jv^ 

3 Nf+l 



/i 

I 1 



C 



2 



Table 1: The critical exponents for mean field theory (MFT) and the Schwinger model 
(SM). Conventions are as in Landau and Lifshitz [5]. 

It is instructive to write down an effective non-local lagrangian from where these values 
of critical exponents can be inferred: 

C = Ar]{A + BT'^y/^fr] + C7]^+Dr]h (2.11) 



Of course, ( p.ll|) is just a shorthand for the values of critical indices obtained and should 
be understood as such. 

In standard theory of second order phase transitions with a nonzero critical tem- 
perature the above 8 critical exponents satisfy 5 universal thermodynamic relations and 
the hyper-scaling relation. The fact that the critical temperature is zero brings about a 
number of modifications compared to the standard theory of phase transitions: 

• The critical exponent /3 refers to the broken phase and therefore cannot be defined. 

• In the strong field limit h ^ t"^^ the partition function depends on the temperature 
as oc exp{ — -^^-^j. Then the critical exponent e becomes singular for t — > and 
cannot be defined. Specifically, for the SM the free energy of the gas of bosons with 
the mass of order /i_ remains exponentially small until t ~ oc m'^ {u = 1 in our 
case). 

• Because Cp = —Td'^^/dT'^ and T = T — Tc (with Tc = 0), an extra power of T 
emerges leading to the Rushbrooke scaling relation 

a + 2/3 + 7 = l (2.12) 



instead of 2. After elimination of (3 and e we have only 3 instead of the usual 5 
relations: 

« + ^ ^ 1 (2.13) 
z/(2-C) = 7- (2.14) 
/i(l + 7-a) = 2v (2.15) 

• The hyperscaling relation follows from the condition that the total free energy is 
of the order of TV/^'^, where V is the spatial volume of dimension d and ^ is the 
correlation length. For Tc = we obtain 

ud=-a (2.16) 

instead of 2 — a. 

The meaning of the hyperscaling relation is that loop corrections to the Green's functions 
estimated from the effective lagrangian ( p.ll| ) are of the same order, as far as powers are 
concerned, as the tree expressions. We emphasize again that one cannot perform serious 
loop calculations with the effective lagrangian (|2.11|) . If trying to do so, logarithmic 



factors would appear which would shift "tree level" values of exponents. 



3 QCD at low temperatures 



The primary interest of the quantity (|1.1|) is that its critical behavior can provide infor- 
mation on the physics of the phase transition in QCD. However, our second remark is 
that the leading infrared behavior of x can be determined exactly in the low temperature 
region. The proper technique to extract it is chiral perturbation theory |l^ . 



Fig. 2. Pseudo-Goldstone loop determining the scalar susceptibility. 
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The leading infrared behaviour is determined by the graph in Fig. 2 involving a loop 
of quasi-massless pseudo-Goldstone bosons. The effective low-energy Lagrangian of QCD 
has the form 

^eff ^ ^F^T,^o,U^)^d^u) + SReTr{A<f/} + ■ ■ ■ , (3.1) 

where U = exp{2i(j)"'t"' / F.^} and 0" are the pseudo-Goldstone fields. The chiral condensate 
is denoted by S = \{qq)o\ (no summing over colors assumed), and Ai is the quark mass 
matrix. The partition function with Ai = diag(m, m, - ■ ■ ,m) has been calculated by 
Gasser and Leutwyler [|l8l . Their expression for the free energy density of lukewarm pion 
gas is 

A^? - 1 A^? - 1 

/ = eo(M.) - -^go{T, M.) + ^^??(T, M.) (3.2) 



where 



g,{T,M^) = -— Jp'dp\n[l-exp{-E/T)], (3.3) 

E = ^/^Tm^, and 

(TM\ 1 f P'dp dgo{T,M^) 
9.{T, M.) = _ = (3.4) 

Here, M^^ stands for the common mass of the N'j — 1 pseudo-Goldstone modes given by 
the Gellmann-Oakes-Renner relation F^M^ = 2mE. 
Substituting in Eq. (P^) the expansion |T^, |2D| 

^°-45^ ~^2~ + "6^' ^^-^^ 

we obtain 

J / ^ 9 TT TT -ir, 



90^ ^ ' 2 "\ 12Nf F2 288N]F^^ 



iV? - 1 „ / 1 T2 
-TM3 1 + 



/V^ — 1 / A2 \ 
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Consider first the case of zero temperature. For tlie infrared singular contribution to 
tlie scalar susceptibility we obtain 

X-^^f^Vlog^ (3 7) 



which was first obtained in |21| from an analysis of the spectral function in the scalar 
channel and the PCAC hypothesis. 

The infrared singular connected contribution to the susceptibility was calculated in |^ 
(One should just substitute 1 for Tr{rt''} = 5"V2 in Eq.(2.13) of Ref.§.). The result is 

Using ( |1.^ ) we immediately obtain the disconnected contribution 

Instead of using the partition function of Gasser and Leutwyler we can calculate the 
susceptibility also using the same technique as in f^. Choosing Ai = diag(m, m, ■ ■ ■ ,m), 
one immediately obtains the vertex 

oy 

mT.Wf\<f>y') = -E^S''\ (3.10) 



which, by evaluation of the diagram in Fig. 2, reproduces the result ( |3.7| ). 

The infrared singular contribution to the susceptibility af finite temperature can be 
obtained directly from (|3.6|) as well 

where the Gellmann-Oakes-Renner relation has been used. The connected contribution 
to the susceptibility was not calculated in P]. However, the zero temperature result of 
Ref. P] can be extended immediately to finite T by making the substitution 

where the sum is over all Matsubara frequencies (po = 27rn//3). The infrared singular part 
comes only from term n = and the result is 

MNf V2^ \FI) 
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This implies that at finite temperature the spectrum of the Dirac operator is nonanalytic 
for small eigenvalues. The disconnected contribution is 

V \ 3/2 



^ JRdis 
AT 



Nj + 2 T 



(3.14) 



The relations (|3.11| , |3.13| , |3.14]) are quite analogous to the relations for the magnetic 
susceptibility for ferromagnets known for a long time. They are also determined by a 
loop of pseudo-Goldstone particles (the magnons) depicted in Fig. 2 and have the same 
behaviorg [0] 



r 



dM 



dH 



T 



H=0 



H 



(3.15) 



a b 
Fig. 3. Two-loop contributions to the scalar susceptibility: a) mass renormalization; b) 
vertex renormalization. 

The relations ( p.ll| , |3.13| , p. 141) hold in the low-temperature region 



(3.16) 



This makes a direct comparison with recent lattice calculations impossible. At not 
so low temperatures, besides the graph in Fig. 2, also higher order graphs in chiral 
perturbation theory contribute. The relevant two-loop graphs are depicted in Fig. 3, but 
the results for the susceptibility can be obtained directly from the partition function as 
well 



f 



T 



'2m 



:*)■"(■ 



1 T^' 



(3.17) 



^Recently, a full nonlinear effective Lagrangian has been constructed [g2| in a way which makes the 
analogy between the theory of ferromagnets and CPT most transparent. 
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The two loop temperature dependence can be absorbed in the one-loop temperature 
modification of the condensate and the pion decay constant as obtained in [|l^ 

FAT) = F,{0) (l - !^'^y (3.19) 

It is not clear what happens to next order in T. Temperature corrections to the 
condensate have been found to three-loop level in pO[. Two-loop effects in M^(T) have 



been extracted in [24, 25 1, but F^{T) at the two-loop level is currently unknown. It is not 
clear whether one can just substitute the temperature dependent S and in ( |3.11| ) to 
all orders in the temperature expansion. The same question can be asked concerning the 
Gellmann-Oakes-Renner relation. Taking into account corrections of order T'^/F^, it also 
holds at nonzero temperatures. Whether it holds at higher orders is an open question 
i2l. 



4 The scalar susceptibility in finite volumes 



The results (|3.7| , |3.11| , p.l3| , p.l4| ) are valid in the thermodynamic limit which means that 



the spatial length of the box where the theory is defined is much larger than the pion 
Compton wave length. However, realistic boxes used in lattice calculations can never 
be made so large if the pion mass is small enough (which is in turn necessary for the 
quantitative analytic predictions to be possible). In this section we consider the opposite 
limit 

A-i < L < ~ (4.1) 

where the susceptibility can be found from the exact results for the finite- volume partition 
function of In this range the susceptibility is expected to be determined by finite 
volume effects and be of order 
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which is much larger than its thermodynamic hmit ~ A^. If oq 7^ 0, the result becomes 
independent of the volume suggesting that it holds for V ^ 00 outside of the range ([4.1|) 
provided that we are sufficiently close to the chiral limit so that S/m ^ A^. This is 
indeed what happens for the pion susceptibility but not for the scalar susceptibility (see 
below). In general, the results of this section have the same status as in Ref.[^. They 
say little about the properties of the theory in the physical infinite volume limit but can 
be used to test the validity of numerical calculations in QCD which are performed for 
finite volumes. 

The theoretically simplest susceptibility is the pseudoscalar (or pion) susceptibility. 
Using the Banks-Casher relation [Q, it can be related immediately to the chiral conden- 



sate, (see eq. ( |1.12| )) 



r = (4.3) 
m 

in the limit that T,mV ^ 1. The validity of this result extends into the thermodynamic 
limit outside of the range ( |4.1| ). 

As follows immediately from ( |l.ll| , prTUD x™"^ can be calculated from average spectral 
density 

p{\) = {u:{\A)) = {^Y.^{\-\^)) (4.4) 

' n 

On the other hand, x*^'^ can be expressed into the connected two-point level correlation 
function 

p,(A, A') = V - p(A)p(A')] . (4.5) 



If, apart from the pairing ±Afc, the eigenvalues of the Dirac operator are uncorrelated 
i.e. if < /(A„)(7(Am) >=< /(A„) >< g{\m) > for n 7^ m and A„, A^ > 0, we have in the 
limit of large V and for positive A, A' 

p,(A,A') = p(A)5(A-A'). (4.6) 

Using the Ua{^) symmetry of the Dirac spectrum the disconnected susceptibility ( |1.10| ) 
can be written as 



r r/A r .y^a^Mhl^^ (A 7) 

Jo "^Vo (A2 + m2)(A'2 + m2) ^^"^^ 



(A2 + m2)(A'2 + m2) 
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which after insertion of the correlation function leads to 



^IRdis ^ MO) ^ S ^ 

m m 



where the second equality is the Banks-Casher relation |^ . 

Because the diagonalization of the Dirac operator induces correlations between the 
eigenvalues, we expect quite a different prediction from the finite volume partition func- 
tion. 

For Nf = 2 the finite volume partition function is known For ^ = 0, 

^ -/i(^S(m„ + m,)) (4.9) 



where V = L^/T in the 4D Euclidean volume. The susceptibilities x'^°'^ ^ind x'^^^ can 
be disentangled by differentiating with respect to different quark masses. Because the 
partition function depends on the quark masses via the sum m„ + rrid we find 

X™"" = 0. (4.10) 

In the limit k, = VHm ^ 1 the disconnected contribution simplifies 

X''' = (4.11) 

In present day lattices the zero mode states are mixed with the much larger number 
of nonzero mode states. Therefore, the partition function is effectively calculated in the 



sector = 0. Existing numerical calculations in the instanton liquid model pO| are also 
done for u = 0. In a sector with fixed topological charge the finite volume partition 
function has been calculated analytically for an arbitrary number of flavors with equal 
mass 



27p. To leading order in k ^, the partition function in the sector of zero topological 



charge is given by 



For K ^ 1 we find 



Z,=o j^jj^. (4.12) 



NfX- + Ny- = ^^. (4.13) 
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In general we have not been able to calculate the connected and disconnected contributions 
to the susceptibility separately. However, for Nf = 2 the partition function is known for 
different quark masses p^. For z/ = we find 



^~ 27r 7o VS(m^ + + 2mumrf cos6')^/2 ' 

which allows us to calculate the connected and disconnected pieces of the susceptibility 
separately. The result for k ^ 1 is = l/2m^\^ and x'^'^ = l/4m^\^. For N j = the 
two contributions to the susceptibility can be obtained from the spectral density and the 
two level spectral correlation function which can be derived from chiral random matrix 
theory. The result for k ^ 1 is x'^"^ = and x'^'^ = l/4m^V". Our conjecture for 
arbitrary Nf and z/ = consistent with all the above results is 

X- = ^ (4.15) 
X''' = ^P'^'i'^)^ (4-16) 



in agreement with the simplest possible flavor dependence consistent with ( [4.131 ). Both 
^con^^^ and F'^^^{k,) approach 1 for k ^ 1, but will in general depend on Nf for finite 
K. The disconnected contribution is suppressed by a factor 1/mVTi with respect to the 
result for uncorrelated eigenvalues. The coefficient oq in ( [4. 21) turns out to be zero which 
is in agreement with the fact that there are no massless scalar particles. 

The quark mass dependence of the scalar susceptibihty has been calculated by lattice 
QCD simulations only for relatively large quark masses 0. In this work with m^V ~ 1 
(in units of the lattice spacing) a quark mass dependence of ~ 1/m is not only found 
at the critical point but also at lower temperatures where the susceptibility levels off at 
a significantly lower value. This result is in between the finite volume prediction ( |4.13|) 
and the result from chiral perturbation theory ( |3.11 ) with mass dependence of 1/m^ and 



l/y/m, respectively. The lattice results for x'^'^ agree with the mass dependence ( |4.8| ) 
suggesting that the eigenvalues of the Dirac operator are only weakly correlated. Two 
different types of correlations can be considered, namely correlations between eigenvalues 
corresponding to different gauge field configurations and correlations between eigenvalues 



obtained from the same lattice gauge configuration. It has been shown that lattice 



eigenvalues show strong spectral correlations (the latter type), but this does not exclude 
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the possibility that correlations between eigenvalues corresponding to different members 
of the ensemble are absent. In random matrix theory it has been shown |^3|] that spectral 
averages and ensemble averages are the same. The exciting possibility that this type 
of generalized ergodicity does not hold for the lattice Dirac eigenvalues deserves further 
attention. 
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